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Abstract The fine structure and the variations of the 
optical pulse shape and phase of the Crab pulsar are 
studied on various time scales. The observations have 
been carried out on 4-m William Hershel and 6-m BTA 
telescopes with APD photon counter, photomultiplier 
based 4-channel photometer and PSD based panoramic 
spectrophotopolarimeter with 1/xs time resolution in 1994, 
1999, 2003 and 2005-2006 years. The upper Hmit on the 
pulsar precession on Dec 2, 1999 is placed in the 10 s - 2 
hours time range. The evidence of a varying from set to 
set fine structure of the main pulse is found in the 1999 
and 2003 years data. No such fine structure is detected 
in the integral pulse shape of 1994, 1999 and 2003 years. 

The drastic change of the pulse shape in the 2005- 
2006 years set is detected along with the pulse shape 
variability and quasi-periodic phase shifts. 
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1 Introduction 

Over the last 30 years the Crab pulsar has been exten- 
sively studied. The reasons for it are clear - it is the 
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brightest pulsar seen in optics, it is nearby and young. 
However, the post popular groups contemporary theo- 
ries of the Crab high-energy emission, the " polar cap" 
CDaug hertv fc Hardind 199611 and "outer gap" ijCheng et al.1 
2000) ones, can't explain the whole set of observational 
data. 

One of the main properties of the Crab emission is 
the very high stability of its optical pulse shape despite 
the secular de crease of the l uminosity, related t o the spin 
rate decrease ()Pacinilll97lt iNasuti et al.lll996j) . 

At the same time the pulsars in general and the Crab 
itself are unstable. The instabilities manifest itself as a 
glitches, likely related to the changes of the neutron star 
crust, timing noise, powered by the collective processes 
in the superfluid internal parts of it, magnetospheric in- 
stabilities, results of the wisps around the pulsar, pre- 
cession, et al. All these factors may influence the optical 
pulse structure and change it on various time scales, both 
in periodic and stochastic way. 

However, it has been found early that the variations 
of the Crab optical light curve, in contrast wit h the radio 
ones, are governed by the Poissonian statistics l|Kristian et al.l 
'1970). A number of observations show the absence of 
non-stationary effects in the structure, intensity and the 
duration of the Crab optical pulses, and the restrictions 
on the regular and stochastic fine structure of its pulse 
on the time scales fr om 3/iS to 500/xs (Beskin ct al. 198^ 
iPercival et alJll993 ), the fluctuations of the pulse inten- 
sitv l)Kristian et aljll970(l . 



Along with the increase of the observational time 
spans and the accuracy of measurements the small changes 
of the optical pulse intensity, synchro nous with the gian t 
radio pulses, have been detected (f Shearer et al.ll20{)3j) . 
Also, the evidence for the short time scale precession of 
the pulsar has been detected by studying its optical light 
curve LCadez et al,.2001 '). 

All this raises the importance of the monitoring of 
the Crab optical emission with high time resolution. 
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Table 1 Log of observations 



Date 


Telescope 


Instrument 


Duration, sec 


Spectral range 


Dec 7, 1994 


BTA, Russia 


Four-color photometer 


2400 


U-fB-l-V+R 






with photomultiplier 






Dec 2, 1999 


WHT, Canary Islands 


Avalanche photo-diode 


6600 


R 


Nov 15, 2003 


BTA, Russia 


Avalanche photo-diode 


1800 


R 


Dec 29, 2005 - 


BTA, Russia 


Panoramic spectro-polarimeter 


48000 


4000 - 7000 A 


Jan 3, 2006 




with position-sensitive detector 







2 Observations 



We analyzed the sample of observational data obtained 
by our group over the time span of 12 years on differ- 
ent telescopes. The details of observations are summa- 
rized in Table Q The equipment used were four-color 
standard photometer with diaphragms based on photo- 
multipliers, fast pho tometer with avalanche photo-diodes 
ijShearer et alJl2003fl and panoramic s pectro-polarimeter 
based on position-sensitive detector ijPebur et al J 12003 
[Plokhotnichcnko et al. 2003). All devices provide the 1/xs 
time resolution. 

For each data set the list of photon arrival times has 
been formed. They have been processed in the same way 
by using the same software to exclude the systematic 
differences due to data analysis inconsistencies. Photon 
arrival times have been corrected to the barycenter of 
the Solar System using the adapted version of axBary 
code by Arnold Rots. The accuracy of this co de has been 
tested with detailed examples provided by iLvne etlil] 
l)2005() and is found to be better than 2/is. 

The barycentered photon lists then have been folded 
using both Jodrell-Bank radio ephemerides ? and our 
own fast-folding based method of timing model fitting. 

The accuracy of timing model is proved to be better 
than at least several microseconds (see Figure QJ, which 
permits to fold the light curve with 5000 bin (6.6 /is) 
resolution. 



3 Phase stability 

We performed the search for timing model residuals us- 
ing two longest continuous data sets of 1999 and 2005- 
2006 years. The data has been divided into the num- 
ber of subsets of fixed length and they have been folded 
separately using the same base epoch. Then the sample 
light curves have been cross-correlated with the standard 
one (which has been derived for each set separately by 
folding the whole data) and its phase shift have been 
derived by fitting the maximum of the cross-correlation 
function with the Gaussian. The results for 1999 year set 
are shown in Figure ^ No evidence for significant devia- 
tions from zero is seen, the phase is consistent with the 
Gaussian noise with 4.1/xs rms in the 10 s - 2 hr time 
range. 



3000 4000 
Time, sec 



Fig. 1 Timing residuals of the Crab pulsar after applying 
second-order timing model (up to second frequency deriva- 
tive). It corresponds to the Gaussian noise with 4.1 fis rms. 



- JD 2453733.5, days 



Fig. 2 Timing residuals of the Crab pulsar after apply- 
ing second-order timing model. The quasi-periodic behaviour 
with characteristic time scale of 0.7 days is seen. 



The data of the last set of 2005-2006 years, however, 
show the significant quasi-periodic variations with ~ 2.5- 
10^'^P rms amplitude. The characteristic time scale of 
the variations is estimated to be roughly 0.7 d. 



Pulse stability of the Crab pulsar 
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Fig. 3 Phased light curves of the Crab pulsar for all obser- Fig. 5 The main pulse peak of the sum of light curves of 
vations sets, scaled to the same pulse height. 1994, 1999 and 2003 years data. 




0.44 0.46 




0.35 
Phase 



Fig. 4 Phased light curves of the Crab pulsar for the three ^ig. 6 The comparison of the peaks of 1999 and 2003 years. 



nights of the Dec 2005 - Jan 2006 set. 



The peaks are shifted vertically for 0.03 for clearance. 



4 Pulse shape 

Due to the presence of significant residuals relative to the 
timing model the pulse profile during the observations 
of 2005-2006 years can't be derived by folding the whole 
data set directly. Instead, we divided the data set into the 
one-hour segments and folded them separately applying 
the time shift corrections to compensate the phase resid- 
uals. The intrinsic phase shift inside each block is less 
than 2 • 10^'*, so the folding with 5000 bins is possible. 
The folded light curves have been co-added. 

All the other data has been folded directly and shifted 
in phase to the same pulse position for the ease of com- 
parison. All pulse profiles are shown in Figure |31 with 
the off-pulse emission subtracted and pulse height scaled 
to the same value. 

The profiles of 1994, 1999 and 2003 years are in a 
perfect agreement with each other. The profile of 2005- 
2006 years, however, deviates from them significantly - 



the pulse remains of the same FWHM while its skewness 
is much smaller, and its shape is nearly symmetric. 

We folded the data of this set for each of three obser- 
vational nights separately using the same method. These 
profiles are shown in Figure ^ There is the significant 
variation of its shape from night to night. Unfortunately 
the low amount of data available do not permit to track 
the profile shape change inside each night and check 
whether it is smooth or whether the shape is correlated 
with the timing residuals. 



5 Pulse fine structure 

For the first three observational set where the pulse pro- 
file is stable we performed the search for the fine struc- 
ture of the main pulse. The data sets has been reduced 
to the same phase base point with precision better than 
half of the phase bin (less than 3.3/is) and the cumula- 
tive light curve has been computed. The peak region of 
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it is shown on Figure No statistically significant de- 
viations from the smooth peak shape is seen. However, 
light curves of 1999 and 2003 years data sets alone (plot- 
ted on FigureEl each show the evidence of fine structure 
on the level of 3-5 sigma (roughly 1 % of the intensity) 
with typical duration of 10-30 /iS. Such details may give 
an evidence of coherent generation of optical emission, 
if the emission generation region is deep enough (deeper 
than 0.1 of light cylinder radius), due to brightness tem- 
perature exceedes 10^^ K. 



other recent change in the neutron star state l)Lvne et alJ 



6 Discussion 

The results of the last data set differs significantly from 
all previous results. We studied carefully the possibility 
of its being the result of some hardware or data process- 
ing problem. 

The data of last set has been acquired using the 
panoramic spectro-polar imeter based on the position- 
sensit ive photon counter l)Debur et al.l2003tlPlokhotnichen' 
12003^ in low-resolution spectral mode. There is no dif- 
ference of the pulse profile in different spectral bands 
(derived using different parts of the spectrum detected). 
The detector behaviour is proved to be linear in the flux 
range used. The data acquisition system, "Quantochron 
4-48" , which records the time of arrival of the detected 
photons, has been checked for short time scale stabil- 
ity by recording the signal from the stationary 100-Hz 
generator. It has been processed and folded in the same 
way as a pulsar one (passing the unnecessary barycen- 
tering step). It shows no distortion of the signal shape 
larger than Ifis. Large scale timing stability of acquisi- 
tion system is ensured by means of 1 Hz and 10 kHz 
frequency signals from GPS receiver. The barycentering 
correction code passed the tests provided in Lvne et al. 
|200 5'1 with the accuracy of 1 /iS. The correctness of the 
radio ephemerides has been checked by performing the 
timing model fitting using our own software, the results 
of phase shift and folding analysis agree with ones based 
on the radio data. There is no small time scale (of order 
of 100 seconds and larger) changes of the pulse profile in- 
side the set with amplitude comparable to the difference 
between the last set and previous ones. 

Taking into account all these arguments we may con- 
clude that the pulse profile change and quasi-periodic 
phase shifts detected in this observational set is most 
likely not related to the hardware or software problems 
of the equipment used. 

The detection of the variations of both the pulse ar- 
rival times and its shape strongly supports the geomet- 
rical interpretation of the effect. It may be described as 
a quasi-periodic change of the pulsar beam orientation 
due to the strong precession commenced suddenly before 
the observations, but after the previous set. It may be 
related to the very strong glitch of the Feb-Mar 2004, or 



7 Conclusions 

We analyzed the data of several sets of optical observa- 
tions with high temporal resolution of the Crab pulsar 
performed by our group over the 12 last years. 

No evidence for short time s cale precession (like 60- 
sec free precession discovered in ICadez" et al.l 1I2OOID ) IS 
detected on the level of 10~^ - 10~^ s~^ pulsar frequency 
variation on 10 s - 2 hours time scale on Dec 2, 1999 (see 
Figure ^ , which corresponds to the precession wobble 
angle to be less than approximately 2 • 10"'^. Also, no 
signatures of short time scale timing noise is seen in this 
data set. 

No significant fine structure is detected in the integral 
pulse profile of 1994, 1999 and 2003 years data set (see 
FigurelU, however, each data set alone show the evidence 
of fine structure on the level of 3-5 sigma, which may be 



a eldLed ito its instability on the time scale of years along 
wlCTi Trie stability of the pulse shape on the same scale. 

We discovered the significant change of the time- 
averaged Crab pulse profile in the Dec 2005 - Jan 2006 
set of observations. The pulse profile also shows the vari- 
ations between the nights. Also, the quasi-periodic phase 
shifts in respect to the second-order timing solution (up 
to second frequency derivative) has been detected in the 
data with amplitude of ~ lOO^s and characteristic time 
scale of 0.7 days. We have not found any hardware or 
software issue able to mimic such pulsar behaviour. These 
results may be interpreted as a geometric effects due to 
the Crab precession suddenly started between our obser- 
vations of 2003 and 2005-2006 years. 
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